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ABSTRACT 

We use frequency-resolved spectroscopy to examine the energy spectra of the 
prominent low frequency Quasi-Periodic Oscillation (QPO) and its harmonic in GX 
339-4. We track the evolution of these spectra as the source makes a transition from 
a bright low/hard to hard intermediate state. 

In the hard/intermediate states, the QPO and time averaged spectra are similar 
and the harmonic is either undetected or similar to the QPO. By contrast, in the softer 
states, the harmonic is dramatically softer than the QPO spectrum and the time aver¬ 
aged spectrum, and the QPO spectrum is dramatically harder than the time averaged 
spectrum. Clearly, the existance of these very different spectral shaped components 
mean that the time-averaged spectra are complex, as also seen by the fact that the 
softer spectra cannot be well described by a disc, Comptonisation and its reflection. 
We use the frequency resolved spectra to better constrain the model components, and 
find that the data are consistent with a time-averaged spectrum which has an addi¬ 
tional low temperature, optically thick Comptonisation component. The harmonic can 
be described by this additional component alone, while the QPO spectrum is similar to 
that of the hard Comptonisation and its reflection. Neither QPO nor harmonic show 
signs of the disc component even when it is strong in the time averaged spectrum. 
This adds to the growing evidence for inhomogeneous Comptonisation in black hole 
binaries. 

While the similarity between the harmonic and QPO spectra in the intermediate 
state can be produced from the angular dependence of Compton scattering in a single 
region, this cannot explain the dramatic differences seen in the soft state. Instead, 
we propose that the soft Compton region is located predominantly above the disc 
while the hard Compton is from the hotter inner flow. Our results therefore point to 
multiple possible mechanisms for producing harmonic features in the power spectrum. 
The dominant mechanism in a given observation is likely a function of both inclination 
angle and inner disc radius. 

Key words: Accretion, accretion discs - X-rays: binaries - X-rays: individual 
(GX 339-4) 


1 INTRODUCTION 

Spectral modeling is a powerful tool in understanding the 
physics of astrophysical sources including accreting black 
holes. By assigning observed/modelled spectral components 
to physical regions of the accretion flow we can learn about 
its geometry and how it changes as a response to changes 
in accretion rate. The most popular physically motivated 
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model for the X-ray spectra of accreting black holes is the 
combination of an accretion disc and an inner hot flow where 
electrons up-scatter the disc photons, giving rise to a comp- 
tonised component. The spectra, especially in the brighter 
states, additionally require the presence of reflection of the 
comptonised photons off some cooler material, presumably 
the accretion disc. Some combination of these components 
can generally fit the observed broadband spectra in all spec¬ 
tral states (for a review, see Done et al.|2007 |. 

However, the data are increasingly showing that more 
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spectral components must be present, especially in the 
brightest low/hard states. The spectra can be matched by 


including extreme inner disc reflection (Wilms et al. 1999 


Fabian & Ross 20101, but alternative models for an addi¬ 


tional component include a warm layer or hot spots on the 
disc ( di Salvo et al.pOOl l, Compton scattering off an addl 
tional thermal or non-thermal electron distribution (Nowak 


et al.|2002 Ibragimov et al.|2005 Yamada et al.|2012|, or a 


component from the radio jet ( Nowak et al.||20TT |. This il¬ 

lustrates the degeneracy in models which are possible using 
only spectral data. 

The degeneracies can be (mostly) broken by using 
additional timing information. Fast timing studies show 
that there are energy and frequency dependent lags across 
the spectrum, both within the Comptonisation compo¬ 
nent and between the disc and Comptonisation component 


( Miyamoto fc Kitamoto||l989 Nowak et al.||1999 Uttley et 
al.|2011| . The lags within the continuum are best explained 
if the Comptonisation is not a single temperature, but is 
softer/cooler at larger radii, and harder/hotter at smaller 
radii. This spatial dependence of the spectrum of the accre¬ 
tion flow, i.e., the inhomogeneous nature of the comptoni¬ 
sation region, is crucial to produce the frequency dependent 
hard lag (soft lead) seen in the data. Longer timescale vari¬ 
ability is generated at larger radii, affecting first the softer 
comptonisation, and then propagates down the accretion 
flow to modulate the fastest variability with a harder spec¬ 
trum produced in the smallest regions close to the black hole 


(Miyamoto &: Kitamoto||1989 Nowak et al.||1999 Kotov et 


al.||2001| [Uttley et al.|200^ 

Since the fast timing results require inhomogeneous 
comptonisation, this means that the broadband continuum 
seen in the time averaged spectrum must be made from 
multiple Compton components. There may be additional 
contributions from the jet and/or reflection and/or inho¬ 
mogeneous disc emission, but the major Compton contin¬ 
uum cannot be a single component. In previous papers we 
have used the technique of frequency-resolved spectroscopy 
( Revnivtsev et ar]|1999 l to track the inhomogeneity of the 
comptonisation region. We used data from the black hole bi¬ 
nary XTE J1550-564 during its rise to outburst as it made a 
transition from the bright low/hard to very high state. We 
found that the fastest variability was always harder than the 
time averaged comptonisation component, consistent with 


this being from the innermost, hardest part of the flow ( Ax¬ 
elsson et al.|2013 (. The strong low-frequency oscillation also 
showed a very similar spectrum to that of the most rapid 
variability, consistent with models where the QPO is from 
Lense-Thirring (relativistic vertical) precession of the inner 
parts of the flow ( Ingram et al.|2009||Axelsson et al.|2014 |. 
These results directly show the inhomogeneity of the comp¬ 
tonisation component. However, the spectral resolution of 
the fast timing modes used for these data was not good, 
with the crucial 8-13 keV data binned into a single point. 

Here, we use better spectral resolution fast timing data 
from a similar transition of the black hole binary GX 339- 
4, and use the frequency resolved spectrum of the QPO and 
its harmonic to constrain the comptonisation components in 
the time averaged spectrum. As before, the QPO spectrum 
is similar to the time averaged comptonisation. However, we 
show for the first time that the spectrum of the harmonic 
can be strongly different from that of the QPO. We suggest 


that there are two major comptonisation regions, a soft one 
where the corona is over the truncated disc, and a hard 
one in the inner regions with no disc. Both are probably 
inhomogeneous. 

The broadband spectra then must be composed of a 
disc, soft Compton (similar to the harmonic) and hard 
Compton continua (similar to the QPO) as well as reflec¬ 
tion. Neglecting the soft compton continua as well as the 
more subtle inhomogeneity in the hard compton continua 
may lead to the more puzzling results from spectral fitting 
such as the very large reflection fraction ( Plant et al.||2014 ( 
and extreme spin (Fabian fc Ross|2010|. 


2 GX 339-4 

The low-mass X-ray binary GX 339-4 was discovered by 
OSO-7 in 1972 ( Marker! et al.|[i973 (, and already this first 
report mentioned strong variability and periods of quies¬ 
cence. It is a transient source, with outbursts occurring on 
average every few years. The compact object is believed 
to have a mass above 7 Mq, and thus to be a black hole 
( Mufioz-Darias et al.|20d8 l. The system inclination has not 
been well determined, but a limit of > 40° is set by assum¬ 
ing the mass of the black hole is less than 20 Mq. Studies 
of the spectral evolution suggest the the inclination is not 
high ( [Mufioz-Darias et al.|[2008| . GX 339-4 is one of the 
most studied black hole transients, with extensive coverage 
by the Rossi X-ray Timing Explorer (RXTE) satellite dur¬ 
ing the 2002, 2004, 2007 and 2010 outbursts. A radio jet has 
been observed from the source ( Gorbel et al.[[2000 (, catego¬ 
rizing it as a microquasar. 

During an outburst, GX 339-4 goes through the canoni¬ 
cal states typically exhibited by black hole binaries. It is also 
the first source in which hysteresis was found in the spec¬ 
tral evolution: the transition from hard to soft state during 
the rising phase occurs at much higher luminosity than the 


soft-to-hard transition during the decline (see, e.g., Nowak 
et al.|2002| ). 

The temporal properties have been studied extensively. 
Following the typical behaviour of black hole binaries, the 
broad-band rapid variability is strong in the hard state and 
suppressed in the softer states. A strong quasi-periodic os¬ 
cillation (QPO) is seen, especially during the hard-to-soft 
transitions, often accompanied by a harmonic. An extensive 
review of the variability features of GX 339-4 can be found 


in Motta et al. (20111. 


3 DATA ANALYSIS 

In this study we use archival data of GX 339-4 from the 
Proportional Gounter Array (PGA; Jahoda et al.|1996 l and 
High-Energy X-ray Timing Explorer (HEXTE: Rothschild 
et al.||19% l instruments onboard the RXTE satellite. While 
there are many RXTE observations made of GX 339-4, only 
a few have both high temporal and spectral resolution, with 
modes suitable for timing often having a single energy bin 
in the 8-13 keV range. We have therefore gone through the 
RXTE archive in search of the best observing modes. 

The selected observations are made in the time period 
2007 February 4-14 (MJD 54135 to 54145), and cover four 
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different stages in the hard to soft evolution. For the PCA 
data we extracted Standard2 spectra for each observation 
applying standard selection criteria. A systematic error of 1 
per cent was added to each bin in the spectra. The energy 
band used in the spectral modelling was 3-30 keV. In the 
case of HEXTE we used Standard Mode spectra from cluster 
B (cluster A was not working at this time), considering the 
energy range 30-100 keV. In the softer states the sensitivity 
does not allow us to go to higher energies. In the harder 
state, we initially attempted to extend the data. However, 
doing so did not change our results and we therefore chose 
to keep the same limit in all observations for consistency. 


3.1 Fourier spectroscopy 


To compare the total spectrum to that of the rapid vari¬ 
ability, we also performed frequency-resolved spectroscopy 
| |Revnivtsev et ^|1999[ |200ip on th e PCA data. Follow¬ 
ing the approach of Axelsson et al. (20131, we extracted 


a light curve for each available channel and constructed a 
power density spectrum (PDS). To obtain the contribution 
from the QPO and harmonic in each energy band, we fit 
them in the PDS using a Lorentzian function and int egrate 
this component (see also Sobolewska & Zycki 20061. This 
gives the relative contribution for each channel, which is fi¬ 
nally combined to make an energy spectrum corresponding 
to each variability component. As in the case of the Stan- 
dard2 spectra, we added a 1 per cent systematic error to the 
frequency resolved spectra. 


4 SPECTRAL MODELING 


4.1 Time averaged spectra 

As described above, the two main components used in ex¬ 
plaining the spectrum of GX 339-4 are the accretion disc 
and a hot inner flow where componisation produces high- 
energy emission. The high-energy photons can then scatter 
off the accretion disc, producing a characteristic Compton 
reflection component. We model the disc component using 
DISKBB, where the main parameter is the maximum temper¬ 
ature Tbb. For the h ard comptonisation component we us e 


the model nthcomp | Zdziarski et al. 


1996 


Zycki et al. 19991, 


parametrized by the asymptotic spectral index F, the seed 
photon temperature To and the electron temperature kT^. 
We set To = Tbb as can be expected if the seed photons are 
the photons from the disc. This assumes that the electrons 
are thermal, though with data extending to higher energies 
is it clear that there are also non-thermal electrons present 


I Joinet et al. 2007 Del Santo et al. 2008 Caballero-Garcfa et 
|al.||2009| ). Howewr, below 100 keV the data are adequately 
modeled by thermal Comptonisation, so we use this descrip¬ 
tion here for simplicity. To model the Compton reflection we 
use XILCONV, a convolution model calculating the reflection 


spectrum using the relxill (Garcfa et al. 20141 code for 


an arbitrary input spectrum. The main parameters are the 
relative amplitude of the reflected component, R, the incli¬ 
nation (which we fix to 50°) and the ionization parameter of 
the reprocessing matter, Finally we include photoelectric 
absorption, keeping the value of the column density frozen 
at Ah = 0.5 x 10^^ cm“^. 


The results of these fits are shown in Table and the 
spectra are plotted in Fig. together with the power spec¬ 
trum in two energy bands. While the fits are acceptable for 
the first two spectra, in the softer states the model clearly 
does not provide a good fit. The residuals indicate that the 
low-energy region is most problematic. 

Overlaid on the spectra are the data points from the 
QPO (blue) and harmonic (green) spectra. They are quite 
similar in the intermediate spectrum, but in the softer spec¬ 
tra they are strikingly different. While the spectrum of the 
QPO appears similar to that of the comptonisation com¬ 
ponent, the harmonic spectrum is dramatically softer. The 
presence of the harmonic in the power spectra means that 
there has to be a component in the time-averaged spectra 
matching its spectrum in these observations. The shape of 
the harmonic spectrum does not match that of the disc and 
we therefore test including an additional soft comptonisa¬ 


tion component. It is modelled using COMPTT (Titarchuk 


19941, with the main parameters being the input photon 


temperature (also here tied to Tbb), the electron temper¬ 
ature fcTaoft and the optical depth r. While COMPTT fur¬ 
ther allows switching between disc and sphere geometries, 
our data do not allow us to discriminate between them, 
and we thus kept the default disc configuration. Though 
the residuals improve, this extra component is not statisti¬ 
cally required for the first two spectra (the significance is 
~ 3 (t); however, there is clear improvement in the two soft¬ 
est spectra (significance > 5cr). Although the soft Compton 
component may be present in all spectra, for our modeling 
we therefore retain it only in the two softest spectra where 
the shape of the harmonic spectrum clearly indicates such a 
component is needed. 

Although the addition of a soft comptonisation compo¬ 
nent improves the fit, the data of the time-averaged spec¬ 
trum cannot constrain all parameters. This is often the case 
when using a complex model, and frequently leads to de¬ 
generacies which are hard to break. In particular the tem¬ 
peratures for the disc and soft Compton components are 
not well determined by the time-averaged spectra alone. As 
the data only start at 3keV the bulk of the disc compo¬ 
nent is not covered and it is therefore not surprising that 
its temperature is difficult to constrain. In the case of the 
soft comptonisation component, it peaks in a region covered 
by the data, but where there are many components present. 
Additionally, there is some inherent degeneracy between the 
electron temperature and optical depth. The uncertainties 
in parameter values are large and some parameters cannot 
be well constrained. This clearly illustrates the difficulty in 
using multiple spectral components. In spite of this, there 
are still some general trends that can be seen in the data. 
For instance, the power-law index evolves from hard to soft 
as can be expected, and the blackbody temperature is higher 
in the softer states. 


4.2 The spectrum of the QPO 

The spectrum of the QPO is somewhat different from to¬ 
tal spectrum, showing less spectral evolution. Studying the 
ratio between the QPO and time-averaged spectrum (Fig¬ 
ure top panel) shows that the QPO spectrum is similar 
to the time-averaged spectrum in the first two observations. 
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Figure 1. Power spectra (upper panels) and fits of our model to the four spectra analysed. The hard-to-soft evolution is clearly seen, as 
is the absence of the harmonic in the first observation. The power spectra in the upper panels are shown for two energy bands, 3—5 keV 
(green) and 10-30keV (blue). In the last two observations, the harmonic is not present in the higher energy band. The lower panels 
show the PC A (cyan) and HEXTE (red) data points overlaid on the total model (solid black line). Also shown are the individual model 
components: disc blackbody (magenta dotted line), main Compton (blue long dashes) and reflection (blue dot-dashed line), as well as 
the data points of the QPO (blue) and harmonic (green) spectra. The bottom panels show the residuals. 


ObsID 

Tbb 

Nbb 

r 

fcTe 

Nc 

R 

log? 

y? DoS 


(keV) 



(keV) 





90235-01-02-02 

n f;Q+0.04 

0+430 

-L- < »_o.02 

36.9l^° 

0 40+0 04 
^•^O-0.07 

0 01+0-13 

o -| 1 +0.09 
'J--‘--L_o.07 

87.0/73 

92035-01-03-00 

'^•'^'^-0.08 

10341^^3 


56tii 

n ^1+014 
^•04_0.10 

0.70l°lg 

q -| q+0.13 

71.7/73 

92035-01-03-03 

n OO+0 05 

1477l«3 


380* 

0 “^o+O-OS 

0-864ri5 

q ^O+0.09 

'J-Uo_0.07 

104/73 

92035-01-03-05 

'^■^^-0.03 

1384l2|2 

0 4Q+0.05 

^•4o_o_o4 

130* 

n 1 q+O-03 

0.72«;« 

q 1 q+0.13 
'^•-^'^-0.09 

128/73 


* Lower limit. 

Table 1. Parameters derived from fitting the time-averaged spectrum (PC A and HEXTE data) with the model 
WABSX (diskbb+nthcomp-I-nthcompxxilconv). Errors indicate 90% confidence intervals. 


but much harder in the last two where the total spectrum 
becomes softer. 

We now test how the spectra of the QPO match our 
model. These spectra have fewer data points and coarser 
energy binning than the average spectra, so it is not possi¬ 
ble to constrain the full model using this data alone. In a 
first step, the best-fitting model for the average spectrum 
is scaled to the level of the QPO spectrum. This provides 
an acceptable fit in the hardest observation, but not for the 
others. In the softest spectra, it is clear that the strong black- 
body and soft Compton component prevent the model from 
fitting the data. 

In order to test for the presence of the disc and soft 
Compton component in the QPO spectra, we allow the nor¬ 
malisation of the blackbody (A^bb), soft Compton component 
(-A^soft), spectral index (P) and amount of reflection {R) to 
vary in the QPO spectrum. Neither the blackbody nor soft 


Compton component is statistically required in the QPO 
spectra, even in the softest states when they are very strong 
in the average spectrum. 

Although the fits converge, we find that even our sim¬ 
plified model is too complex to be uniquely constrained by 
the data. The amount of reflection appears to be consis¬ 
tently high in the QPO spectra, with the maximum allowed 
reflected fraction of 1.0 within the error range. The spectral 
index is similar to that of the average spectrum even in the 
softer states, which seems in contradiction with Fig. A 
closer look reveals that this is coupled to the strong reflec¬ 
tion. The reflection “hump” at higher energies is balanced 
by a steep spectral index to give a stronger comptonisa- 
tion contribution at lower energies. The result is a fairly fiat 
overall spectrum, as seen in the lower panels of Fig We 
stress that a weak soft Compton component would also give 
a contribution at low energies; there is indeed a negative 
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Figure 2. Ratio between the QPO and average spectrum (top), 
QPO and main comptonisation with reflection (middle panel) and 
harmonic to QPO (bottom). All ratios are normalised to 1 at 5 
keV. The colours show the evolution during the outburst. From 
hard to soft: cyan, red, blue, green. 

correlation between the strength of this component and the 
spectral index. 

While we cannot rule out a weak contribution from the 
disc or soft Compton components, the simplest model which 
can match the QPO spectrum is the combination of main 
comptonisation and its reflection. We find that these two 
components alone are able to give good fits to all four ob¬ 
servations. This is also evident in the middle panel of Fig.[^ 
showing the ratio of the QPO spectrum to the main comp¬ 
tonisation and reflection from the time-averaged spectrum. 

4.3 The spectrum of the harmonic 

The spectra of the harmonic show dramatic differences be¬ 
tween the observations. In the low/hard state (ObsID 92035- 
01-02-02), no harmonic is present, but the spectra in the 
other states are shown in the lower panels of Fig. The 
three observations are all from different stages of the hard in¬ 
termediate state, yet the harmonic spectrum varies greatly. 
In the hardest of the three (ObsID 92035-01-03-00) the har¬ 
monic is stronger than the QPO, yet shows a very similar 
spectrum. In the two softer observations, the harmonic is 
weaker and its spectrum very different from either the av¬ 
erage spectrum or the QPO. The greatest change is seen at 
high energies, with the harmonic disappearing in the power 




10 100 10 100 
E [keV] E [keV] 


Figure 3. Spectral fits for the two softest observations when in¬ 
cluding a soft comptonisation component in the model. The upper 
panels show the PCA (cyan) and HEXTE (red) data points over¬ 
laid on the total model (solid black line). Also shown are the indi¬ 
vidual model components: disc blackbody (magenta dotted line), 
soft Compton (green short dashes), main Compton (blue long 
dashes) and reflection (blue dot-dashed line). The data points of 
the harmonic (green) spectra are overlaid on the scaled-down soft 
Compton component (lower solid black line). The bottom panels 
show the residuals. 

spectra above ~ lOkeV. At first glance, the harmonic spec¬ 
trum in the two last observations appears similar to a black- 
body spectrum, although the peak is at higher energy than 
the disc. We therefore attempt to fit them using a single¬ 
temperature blackbody, e.g., corresponding to the harmonic 
arising in a confined region at the inner boundary of the disc. 
However, this does not provide a good fit - the spectrum of 
the harmonic is too broad. 

As the presence of the harmonic spectrum motivated us 
to include the soft Compton component in the time-averaged 
spectra of the softer states, it may provide a means to con¬ 
strain this component. To test this, the harmonic spectrum 
is fit together with the averaged one (PCA-I-HEXTE). The 
model used is the same as before, but the harmonic spec¬ 
trum is fit using the soft comptonisation component only 
(i.e., wabsxcomptt). The parameters of this component 
are tied to be the same in the average spectrum, only al¬ 
lowing the normalization to differ. The results are shown in 
Figure and the parameters given in Table 

The fit parameters are in general compatible with those 
of the average fits alone, but constraining the soft Comp¬ 
ton component with the harmonic allows more robust fits 
with smaller uncertainties. However, the degeneracy be¬ 
tween electron temperature fcTsoft and optical depth r in 
the soft comptonisation prevent us from constraining both 
parameters. We therefore freeze the temperature at its best¬ 
fitting value when performing the error estimations. Having 
done this, the added data from the harmonic allow us to 
constrain the remaining parameters, which was not possible 
using time-averaged data alone. 


5 DISCUSSION 

Comparing the power spectra of the two energy bands in 
Fig.[T] in the hard state they are almost identical. The vari¬ 
ability here is thus more or less independent of energy, and 
likely we are probing the same component (i.e., the main 
comptonisation) in both bands. A completely different pic- 
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ObsID Tbb A^bb kT^oit r N^aft T fcTe Nq R log$ X^Dof 

(keV) (keV) (keV) 


92035-01-03-03 O-OOlo;®* 2400124^0 l-lSloJg 0.27to:?6 2.191°;“ 39' 

92035-01-03-05 0.68lg;°® 4060l°3° 0.78l°-°g 0.26lg;°° 2.191°;“ 36' 


0.44lg;“ 0.75lo;ig 3.05lb;3g 74.0/74 
0-20lp°4 0.98l°;°g 2.99lg;°'[ 61.9/75 


^Parameter frozen at this value during error calculations. 
' Lower limit. 


Table 2. Parameters derived from fitting the time-averaged spectrum (PCA and HEXTE data) simultaneously with the spectrum of 
the harmonic (except for ObsID 92035-01-02-02). The spectrum of the harmonic is fit using only the soft comptonisation. Errors indicate 
90% confidence intervals. Note that the reflection fraction was constrained to ij < 1 in the fits. 


ture is instead seen in the softer observations. In these cases, 
the shapes of the PDS match up until ~ 1 — 2 Hz, and then 
start to differ. This pattern cannot be explained by a sin¬ 
gle varying component. If the strength of this component is 
different between the bands, this would led to a dilution of 
the variability, and thus a shift in the PDS, but the over¬ 
all shape would be the same. Neither can the behaviour be 
explained by “pivoting” of the spectral index, as this would 
also produce a shift in the relative strength (depending on 
the location of the pivot point) without altering the shape of 
the PDS. Instead, the observed behaviour points to an addi¬ 
tional component with a characteristic variability timescale 
of ~ 2 Hz. This is too rapid for the cool disc, and indicates 
the presence of a separate hot/geometrically thick compo¬ 
nent in addition to the main comptonisation. It is thus clear 
already from these considerations that the comptonisation 
region is inhomogeneous. 

Our result that the QPO spectrum is similar to 
the main Comptonisation agrees with previous results 
from XTE J1550-564, though in GX 339-4 the QPO is 
subtly softer than the main Comptonisation whereas in 
XTE J1550-564 it is subtly harder. However, the behaviour 
of the harmonic is very different between the two objects. 
In J1550, the harmonic was always softer than the QPO, 
but it could be fit with a comptonisation spectrum with the 
same electron temperature as for the QPO and main Comp- 


lack of wind features: Ponti et al. (2012l). Thus we suggest 


tonisation component (Axelsson et al. 20141. By contrast. 


in GX 339-4, the harmonic seen in the softer spectra clearly 
requires a much lower electron temperature than the main 
Comptonisation component. 

Precession of the hot inner flow can produce both a 
QPO and harmonic (from the angle dependence of Comp¬ 
ton scattering and the non-spherical nature of Compton re¬ 
gion) so this predicts that these have spectra which are sim¬ 


ilar to that of the Comptonisation component ( Ingram et 
al.||2009 Veledina et al.|[2013L Subtle differences in spectra 


can be produced by the Compton region being inhomoge¬ 
neous, and the QPO and harmonic being weighted to dif¬ 


ferent radii (Axelsson et al. 20141. However, the dramatic 


difference of the harmonic in the softest states of GX 339- 
4 seems too marked to be incorporated into this picture. 
Instead, it suggests that there is a another mechanism to 
produce the harmonic, and that this mechanism is more ev¬ 
ident in GX 339-4 than in J1550. Yet the two systems are 


thought to be very similar in binary parameters (Munoz- 


Darias et al. 20081, with the only difference being inclina¬ 


tion angle. J1550 is known to be at a rather high inclination 
angle (74.7±3.7 degrees; Orosz et al. pMTt , while GX 339-4 
is viewed at intermediate inclination (as inferred from both 
its low disc temperature: Munoz-Darias et al. (20131, and 


that the visibility of the additional harmonic signal depends 
on inclination, being more evident at low inclination angles. 
This in turn suggests that it connects to the disc, as is also 
supported by it being seen in the time averaged spectrum 
mostly in states where there is a strong disc component. 


Motta et al. (20151 find differences between the strength 


of QPOs as a function of inclination, with different depen- 
dance for different types of QPOs. However, in our study 
we see no big differences in the behaviour of the QPO be¬ 


tween GX 339-4 and XTE J1550-564 (studied in Axelsson 


et al. 20141. The big difference is in the harmonic, which 


in the harder states of GX 339-4 is similar to J1550, but in 
the softer states changes drastically. This means that incli¬ 
nation alone cannot be the cause for the difference in the 
harmonic spectrum. Studies comparing the coherence as a 
function of frequency for type G QPOs in several black hole 
sources ( Rao et al.|2010 Pawar et al.|2015 \ have found that 
there are differences between the different harmonics. Some, 
like the fundamental, are frequency modulated while others 
appear amplitude modulated. This could tie into the results 
found here, indicating more than one different mechanism 
behind the harmonics. 

The similarity found between the harmonic spectrum 
and the soft comptonisation component in the average spec¬ 
trum of GX 339-4 clearly suggests they are related. Yet it is 
not obvious how such a component could produce the har¬ 
monic, or why the other emission components do not con¬ 
tribute. If the soft Gompton component is the outer edge of 
the inner flow (the region closest to the disc), it would be 
expected to overlap the disc during the softer states. This 
prevents precession, and thus we would not expect the re¬ 
gion to be present in the harmonic (or QPO) spectra at 
all. As noted above, the QPO spectra do not require a soft 
Compton component, suggesting that this region is indeed 
located further out in the flow. 

Another picture could be that the harmonic instead 
arises because of the disc. If the variability is a vertical mode 
in the inner flow, the stabilizing disc could in practice cut the 
wavelength in half and thereby cause the frequency to dou¬ 
ble. This would occur preferentially in the outer edge of the 
hot flow, matching the geometry inferred by, e.g.,|Yamada et| 


al. (2012 I. The QPO picks out the flow which can vertically 
precess i.e., the hot inner flow. As this comes past the disc 
it drives a wave into the softer comptonisation region which 
bounces off the midplane disc, so the harmonic picks out the 
spectrum of the softer comptonisation region. A doubling of 
one of the characteristic frequencies in Cyg X-1 was also 
reported by Axelsson et al. (20061 during the state tran¬ 


sition, suggesting this could be a common behaviour when 
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the hot inner flow and accretion disc overlap. Although emis¬ 
sion from the soft Compton component should be reflected 
on the disc, the resulting feature would be shifted down in 
energy and not appear in the RXTE bandpass, explaining 
the match between harmonic and direct emission. 

The need for a soft Compton component has also been 
found in other sources, such as Cyg X-1. In a systematic 
study of spectral evolution in the hard and intermediate 
states, Ibragimov et al. (20051 found that the optical depth 


of the soft Compton components decreases for the softer 
states. However, Makishima et al. (20081 find no significant 
difference in optical depth of their soft Compton component 
between the hard and hard intermediate state. Although the 
optical depth is lowest in our softest observation, since we 
only clearly see the component in two spectra we cannot 
make any definite statements. 

The geometrical origin of the harmonic, as well as an 
origin in the inner disc or overlap region, is strongly depen¬ 
dent on viewing angle and thereby inclination. This could 
explain the differences seen between GX 339-4 (low to mod¬ 
erate inclination) and J1550 (high inclination). In addition, 
both mechanisms also depend on the inner radius of the 
accretion disc. Nevertheless, the model proposed here pre¬ 
dicts that the geometric origin of the harmonic should be 
strongest in high-inclination sources, whereas the mecha¬ 
nism behind the softer harmonic spectra (presumably re¬ 
lated with the inner disc region) is seen for low and moderate 
inclination sources. This can be tested with future system¬ 
atic studies. 

The results also demonstrate that the geometry is likely 
complex. Comparison of the QPO and average spectra both 
here and in Axelsson et al. (20141 show that the hard Comp¬ 


ton region is inhomogeneous (see also Hjalmarsdotter et al. 


2015). It is therefore not difficult to imagine that the same 


would apply also to the soft Compton region. Our simpli¬ 
fied models are not sufficient to capture these details, and 
the available data are not able to constrain more complex 
combinations of components. We have assumed a single tem¬ 
perature of both the input seed photons as well as the comp- 
tonising electrons; in reality, there will be a distribution of 
temperatures in both cases. In an overlap geometry the disc 
temperature observed directly is also likely to be lower than 
that seen by the comptonising region covering the inner re¬ 
gions of the disc. Despite these difficulties, we find that using 
our model we can build a coherent picture to explain the av¬ 
erage spectra together with those of the QPO and harmonic, 
as well as their evolution. This is an important first step to¬ 
wards a physical picture of both spectra and variability in 
black hole binaries. 


6 SUMMARY AND CONCLUSIONS 

We have studied the spectrum of the QPO and harmonic in 
GX 339-4 using the RXTE observations with the best com¬ 
bination of spectral and temporal resolution, and compared 
them to the time-averaged spectrum. Our first result, by 
comparing power spectra in the soft and hard bands, is that 
the comptonising region is inhomogeneous. This is clear evi¬ 
dence for an additional low-energy component being present 
in the spectrum, at least in the softer spectra. 

The QPO spectrum is harder than the time-averaged 


one in the softer states and can be well fit using two compo¬ 
nents: thermal comptonisation together with Compton re¬ 
flection off an ionized disc. This matches predictions tying 
its origins to precession of the inner flow. 

The drastic change seen in the spectrum of the har¬ 
monic indicates that more than one mechanism must be able 
to produce it. For harder spectra, the similarity with the 
QPO spectrum suggests the same geometrical origin as in 
J1550. For the softer observations the spectrum of the har¬ 
monic matches the soft Compton component present in the 
time-averaged spectrum, providing a means to constrain this 
component. The absence of the other spectral components 
in the data of the harmonic places its origin, and thereby 
that of the additional low-energy component, at the outer 
edges of the hot flow. We speculate that it is related to the 
overlap region, with the harmonic arising as the disc pen¬ 
etrates the flow. Which of these mechanisms dominate is 
likely a complex function of inclination angle and inner disc 
radius. We propose systematic studies of multiple sources of 
varying inclination as a means to test this. 

Our results clearly show the power of combining tem¬ 
poral and spectral information. By adding the information 
from the variability, the different components in the spec¬ 
trum can be constrained - something which is not possible 
by spectral fitting alone. A natural extension is to combine 
this with phase-resolved techniques (Ingram et al, in prep.), 
giving information also on possible modulation of compo¬ 
nents with orbital phase. 
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